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Resumen 

Esta tesis tiene el proposito de estudiar como se modifican algunos de los observables as- 
trofisicos de Estrellas de Quarks debido a la presencia del campo magnetico. 

Para hacerlo, es necesesario obtener la Ecuacion de Estado (EdE) y estudiar la estabilidad 
de la Materia Extrana (formada por quarks u, d y s) magnetizada densa y frfa en equilibrio 
estelar (equilibrio (3, conservation del numero barionico y neutralidad de carga). Trabajaremos 
utilizando el modelo fenomenologico de Bag del MIT. 

La estabilidad de la Materia Extrana de Quarks Magnetizada (MEQM) se estudia teniendo 
en cuenta la variation de los parametros mas importantes del modelo: la masa del quark s, la 
densidad barionica, el campo magnetico y el parametro de Bag. Los resultados obtenidos se 
comparan con los de la materia normal de quark magnetizada (solo quarks uy den equilibrio j3) 
asi como con los de la Materia Extrana de Quarks (MEQ). Se obtiene que la energfa por barion 
decrece con el aumento del campo magnetico lo que implica que la MEQM es mas estable que 
la MEQ. 

Las EdE previamente obtenidas se utilizan para obtener configuraciones estables de Estre- 
llas Extranas Magnetizadas, comprobando que el campo magnetico contribuye a disminuir la 
relation Masa-Radio (M-R) de la estrella. 



Abstract 

This thesis aims to study how amending certain astrophysical observable of quark stars due 
to presence of magnetic field. 

To do this we need to obtain Equation of State (EOS) and consider the stability of Strange 
Quark Matter (made up of quarks u, d and s) cold dense and magnetized in stellar equilibrium 
(j3 equilibrium, conservation of the baryonic number and charge neutrality). We will work using 
the phenomenological MIT Bag model. 

The stability of the Magnetized Strange Quark Matter (MSQM) is studies taking into 
account the variation of parameters from the model: s quark mass, baryonic density, magnetic 
field and the Bag parameter. Results obtained were compared with those of magnetized normal 
quark matter (only u and d quarks in (3 equilibrium) as well as the Strange Quark Matter 
(SQM). It is found that the energy per baryon decreases with the increasing magnetic field 
which implies that the MSQM is more stable than SQM. 

The Equations of State previously obtained are used to obtain stable configurations of 
magnetized strange stars checking that the magnetic field helps to reduce Mass-Radius (M-R) 
ratio of the star. 
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Introduccion general. 



Los Objetos Compactos (OC) son los remanentes que quedan tras la muerte de una estrella. 
Dependiendo de sus caracteristicas se clasifican en Enanas Blancas (EBs), Estrellas de Neu- 
trones (ENs) y Huecos Negros (HN) pQ. En estos objetos, la materia esta sujeta a condiciones 
extremas (altfsimas densidades, altas temperaturas y campos magneticos muy intensos), ello 
hace que sean considerados laboratorios naturales para probar teorias, que dado el desarrollo 
tecnologico actual, son imposibles de comprobar en laboratorios terrestres [2]. 

En una estrella comun existen dos fuerzas que se contraponen de forma tal que se llega a 
un estado de equilibrio mecanico. Estas fuerzas son: la gravitacional, que tiende a comprimir 
la materia hacia el centro; y la presion ejercida por la radiation desprendida en las reacciones 
de fusion, que tiende a expulsar materia hacia el exterior de la estrella. 

Durante millones de anos la estrella quema el combustible nuclear manteniendo un esta- 
do de equilibrio, ubicandose en la llamada secuencia principal del diagrama de Hertzsprung- 
Russell [3J. Cuando en la estrella cesan los procesos de fusion, la gravedad comprime la materia 
y es la presion del gas degenerado de electrones, en el caso de las EBs; y de neutrones, en el de 
las ENs, lo que garantiza la estabilidad de la estrella. 

La complejidad de Enanas Blancas y ENs requiere que para explicarlas dominemos toda la 
Ffsica que conocemos: la Relatividad General porque en muy poco espacio esta concentrada 
mucha materia y el espacio tiempo se deforma de manera apreciable, la interaction electro- 
magnetica porque hay interaction de particulas cargadas que compensan la carga de la estrella 
resultando una estrella neutra, la interaction electrodebil porque hay constante decaimiento 
de neutrinos y la fuerte porque o hay quarks deconfinados o hay hadrones^ conformados por 
quarks. 

Muchas interrogantes surgen cuando estudiamos este tipo de objetos. En particular para 
las ENs los grandes campos gravitacionales hacen que las densidades en ellas sean del orden de 
la densidad del nucleo atomico y hasta varias veces mayor (la masa del Sol esta concentrada 
en un radio de ~ 10 km) [I]. En la medida que nos acercamos a su centro, pudieran ocurrir 
procesos aun mas exoticos que puedan producir que los quarks "salgan" de sus neutrones, y se 
forme entonces una Estrella de Quarks (EQs). En realidad serfa una Estrella Extrana porque 
aparece ademas de los quarks up(u) y down(d) otro mas masivo: el quarks strange o extrano 
(s), debido a procesos de desintegracion /3 donde se transforman quarks u y d en s [5]. 

Estos objetos serian muy diferentes a las ENs, podrfan ser mas pequenos y compactos, tener 

1 Se dividen en bariones(formados por 3 quarks) y mesones (pares quark-antiquark) . 
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unos 6km de radio y rotar mas rapido. La materia estarfa cohesionada debido a la interaction 
fuerte y no a la gravitatoria como ocurre en las ENs [6]. 

Por tanto, conocer que compone el interior de una ENs o afirmar que la existencia de 
una EQs es una realidad, esta condicionado a responder cual es el estado fundamental de la 
materia a muy altas densidades y con ello responder a la pregunta, desde la fisica de particulas 
elementales: ique le ocurre a la materia cuando es "comprimida" hasta alcanzar densidades por 
encima de la densidad nuclear (no = 0,16 fm~ 3 ) ^ 

Hoy en dfa, la teorfa aceptada para describir el comportamiento de las particulas subatomi- 
cas y sus interacciones es el Modelo Estandar. En este modelo se utiliza el aparato matematico 
de la Teorfa Cuantica de Campos para unificar el Modelo Electrodebil (interaction electromag- 
netica y debil unificada) y la Cromodinamica Cuantica (CDC) en una teorfa que describe las 
interacciones de todas las particulas observadas experimentalmente. 

Segun el Modelo Estandar, la materia esta compuesta por doce fermiones: seis quarks y seis 
leptoneajjy sus antipartfculas), ademas estan los bosones mediadores de las interacciones. En 
la Tablalllse muestran algunas de las propiedades de los quarks [7]. 



Quarks 


Masa (MeV) 


Carga(e = — 1) 


u (up) 


1,5-4,0 


+2/3 


d (down) 


4-8 


-1/3 


s (strange) 


80 - 130 


-1/3 


c (charm) 


1150 - 1350 


+2/3 


t (bottom) 


4100 - 4400 


-1/3 


b (top) 


170900 ± 1800 


+2/3 



Tabla 1: Masas y cargas de los quarks [8]. 

Los quarks son los "ladrillos" que junto con los electrones forma la materia tal y como la 
conocemos: los nucleones (neutrones y protones), que estan formados por quarks u y d y junto 
a los electrones forman atomos y estos las moleculas. El hecho de que ellos esten unidos dentro 
de los nucleones se debe a la interaction fuerte o de color (CDC) que opera a cortas distancia 
~ 10~ 13 cm. 

La CDC describe las interacciones entre los quarks y los gluones, particulas mediadoras de la 
interaction fuerte, analogas a los fotones (mediadores de la interaction electro magnetica) . Entre 
las principales caracterfsticas que se derivan de la CDC estan el confinamiento y la libertad 
asintotica que dan cuenta de que no hayan quarks libres en la Naturaleza, esten confinados a 
los hadrones y se comporten como libres dentro de ellos. 

Sin embargo la CDC concibe que en condiciones extremas: alta temperatura o alta densidad, 
pueden existir quarks libres [HI HP] . Por tanto, pensar que el interior de una ENs pudiera ser el 
escenario real de la existencia de este tipo de materia es totalmente plausible. En el diagrama 



2 lfm = 10" 13 cm 
3 neutrinos y electrones 
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DENSIDAD BARIONICA 



Figura 1: Esquema del diagrama de fase de la CDC 



de fase de la CDC Fig (JTj) se representa la zona en que "habitan" los OC, es decir, bajas 
temperaturas (T ~ 30MeV) y altos potenciales qumiicos (jj, > 300 MeV). 

Los radios y masas de las EQs o las ENs se pueden medir indirectamente. Ellas emiten 
radiacion Gamma y rayos X que se detectan por instrumentos colocados en satelites |T1J. A 
estas estrellas se les asocian ademas campos magneticos muy intensos que se explican por la 
amplification del campo de las estrellas progenitoras despues de la explosion de la Supernova. 
Estos campos son extraordinariamente grandes: de 10 15 Gauss (el campo terrestre es de 0.5 
Gauss) en la superficie, y 10 18 Gauss (G) en el interior |3]. La NASA en el 2002 reporto 3 
fuentes candidatas a Estrellas Extranas |12] . 

En esta tesis nos proponemos estudiar la Materia de Quarks (MQ) en presencia de campos 
magneticos intensos en equilibrio estelar y los observables astrofisicos que se obtienen para EQs 
formadas por MQ magnetizada. 

Este trabajo le da entonces continuidad a [131 EJ US] donde fueron estudiadas las propie- 
dades termodinamicas del gas de quarks magnetizado porque, no solo se estudia la estabilidad 
de la materia de quarks y las EdE en condiciones realistas, sino tambien se investigan las im- 
plicaciones macroscopicas que tendrfa la existencia de EQs formadas por este tipo de materia 
magnetizada. Esto ultimo se hace tomando en cuenta los efectos de la Relatividad General en el 
equilibrio hidrodinamico de la estrella a traves de las ecuaciones Tolman-Oppenheimer-Volkoff 
(TOV) PS]. 

Para ello usaremos el Modelo Fenomenoloeico MIT Bag Model que permite describir 
el confinamiento de los quarks sin necesidad de usar la Cromodinamica Cuantica CDC que 
introduce grandes complicaciones matematicas. En dicho modelo los quarks son considerados 



Massachusetts Institute of Technology. 
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libres dentro de un volumen (Bag) y es precisamente este quien reproduce el confinamiento. El 
equilibrio estelar significa que consideraremos que transcurre suficiente tiempo como para que 
ocurran procesos de desintegracion dentro de la estrella, mientras que la misma se mantiene 
electricamente neutra y conserva el numero de bariones en su interior. 

Un estudio asf tiene importancia no solo para la Astrofisica sino tambien para la Ffsica de 
Partfculas. Buscando concordancia con las observaciones podrfamos restringir parametros de 
los modelos teoricos que influirfan en las propiedades de la ffsica del micromundo, hasta que 
lleguen experimentos terrestres que serian los mas conclusivos. 

En la tesis nos proponemos concretamente los siguientes objetivos: 

1. Obtener la EdE, relation entre la presion y la densidad de energfa de la MEQ y la materia 
normal de quarks (materia electricamente neutra formada por quarks Myden equilibrio 
con electrones), en presencia de campos magneticos intensos en equilibrio estelar. 

2. Estudiar la estabilidad de la MEQM en equilibrio estelar. Probar que la energfa por 
barion para la MEQM es menor que para la MEQ y para la materia normal de quarks 
magnetizada. 

3. Resolver las ecuaciones de TOV usando las EdE previamente obtenidas y asf determinar 
las configuraciones estables de masa-radio para las estrellas de quark extranas magne- 
tizadas. Obtener ademas la masa barionica y el corrimiento al rojo gravitacional para 
EQs. 

La tesis esta organizada de la siguiente forma: Los Capitulos [T] y [2] estan dedicados a 
introducir los aspectos fundamentales referente a los objetos compactos y a las propiedades de 
la materia de quarks respectivamente. En ellos solo pretendemos dar a conocer los principales 
avances en estas areas, por lo cual profundizaremos solo en los aspectos que nos sean necesarios 
en los demas capitulos. En el Capitulo|3]exponemos las principales caracterfsticas de la materia 
de quarks magnetizada, obtenemos las EdE y estudiamos la estabilidad de la materia de quarks 
magnetizada en equilibrio estelar. En el Capitulo [4] obtenemos las configuraciones estables 
de masa y radio para EQs, asf como la masa barionica y el corrimiento al rojo gravitacional. 
Posteriormente se exponen las conclusiones y recomendaciones de la tesis. 



Capitulo 1 

Objetos Compactos: Estrellas de 
Neut rones, Estrellas Ext r anas. 



1.1. Caractensticas generales. 



En el Universo existe una gran variedad de fenomenos que despiertan el interes y el asombro 
de los seres humanos, sin dudas las estrellas constituyen uno de los escenarios mas atractivos 
para la investigation. 

La vida de una estrella pasa por diferentes etapas. Todo comienza cuando las nubes de 
gas estelar comienzan a agruparse gravitacionalmente, la temperatura aumenta y comienzan a 
fusionarse los elementos ligeros como el hidrogeno. En estas reacciones de fusion se van formando 
elementos mas pesados. Cuando el combustible nuclear se agota, la estrella se transforma en 
una Supernova, su destino final depende en gran medida de la masa inicial de la estrella |17] . 



En la Fig (1.1) se resumen los posibles destinos de una estrella en dependencia de su masa, y 



sus principales caractensticas se resumen en la Tabla |1.1 

Entre los posibles destinos finales de las estrellas, nos centraremos en las ENs, debido a que 
como ya mencionamos con anterioridad ellas podrfan tener en su interior una fase formada por 
MEQ formando o una Estrella Extrana o una hfbrida. 



Objeto Compacto 


Masa a (M) 


Radio 6 (R) 


Densidad Media ( g cm 3 ) 


Enana Blanca 


< 1,4M 


~ io- 2 i? 


<io 7 


Estrella de Neutrones 


~ 1 - 3M 


~ icr 5 i? 


<io« 


Hueco Negro 


Arbitraria 


2GM/c 2 


~ M/R 3 



a Masa del SolM 
6 Radio del Solifc 



■0 



1,989 x 10 dd g 
6,9599 x 10 10 



cm 



Tabla 1.1: Caractensticas de los Objetos Compactos [3]. 
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Figura 1.1: Diagrama que muestra el destine- final de las estrellas en dependencia de su masa [17j. 

1.2. Estrellas de Neutrones. 

Las ENs fueron predichas por Landau poco tiempo despues de haberse descubierto el neutron 
en el ano 1932, desde entonces mucho se ha avanzado en el estudio de sus propiedades. Una 
ENs es el remanente de una estrella supergigante roja que al agotarse los elementos necesarios 
para los procesos de fusion comienza a colapsar, eyectando las capas externas en una explosion 
conocida como una Supernova tipo II. La materia que queda en el centro de la estrella se 
encuentra sometida a campos gravitacionales muy elevados. Debido a las altas densidades, los 
electrones son capturados por los protones formandose neutrones en una reaction /3 inversa 
(p + e~ — > n + v e ). La presion del gas degenerado de neutrones es la que frena el colapso de la 
estrella equilibrando de esta manera la fuerza gravitational pi]. 

Las ENs estan formadas por capas en las cuales podemos encontrar distintas familias de 
partfculas. Se pueden identificar cinco capas fundamentales pi]: el nucleo interior y exterior, la 
corteza, la envoltura y la atmosfera. En la atmosfera y la envoltura se encuentra una cantidad 
despreciable de masa de la estrella, estas capas estan formadas principalmente por electrones, 
nucleos y atomos. La corteza se extiende por unos 2 km, esta compuesta por nucleos que varian 
de acuerdo a la densidad, van desde el nucleo de 56 Fe para densidades de 10 6 gcm -3 hasta 
nucleos con numero barionico A ~ 200 para densidades n ~ no/3. A medida que la densidad 
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continua aumentando en la corteza hay toda una serie de cambios de fase hasta llegar al nucleo 
donde se concentra el 99 % de la masa. El nucleo exterior es una mezcla de nucleones, electrones 
y muones. En el nucleo interior se espera que exista una transition de fase a materia deconfmada 
de quarks. 

Una caracterfstica importante de las ENs son sus elevados campos magneticos (B ~ 10 12 — 
10 14 G), el mecanismo mediante el cual surgen estos campos no se comprende del todo. El modo 
mas simple de explicar su surgimiento es producto de la compresion del flujo magnetico durante 
el colapso del nucleo de la estrella progenitora, con esto se logra explicar la formation de campos 
entre 10 11 - 10 13 G. 

Hay otro efecto que ocurre dentro de algunas ENs que pudiera explicar campos magneticos 
mil veces mayores. Dentro de ellas el "gas de partfculas"circula por conveccion, es decir, la 
temperatura no es uniforme luego hay partes del gas que estan mas calientes que otras. Como el 
gas es un buen conductor electrico cualquier lmea de campo magnetico lo arrastra moviendolo 
con el. De esta manera el campo puede amplificarse. Este efecto se conoce como el Efecto 
Dinamo [TBI [T9] y se piensa que es el que origina el campo magnetico en todas las estrellas e 
incluso en planetas; esta presente durante toda la vida de las estrellas de una manera "sedada" 
pero para las EN cuyo nucleo rota rapidamente la conveccion llega a ser muy violenta. La 
primera simulation numerica hecha para ENs a partir de un modelo con conveccion puede verse 
en |2Uj y los resultados muestran que si una Estrella Neutronica naciente tiene una temperatura 
de 3 x 10 9 K el ffuido caliente circula en diez milisegundos y luego cesa. Esta furiosa conveccion 
explica el magnetismo de este tipo de estrella de neutrones [18]. En otras palabras, este modelo 
nos explicarfa que si una estrella naciente de neutrones entrega el 10 % de su energfa cinetica 
al campo magnetico este puede llegar a 10 15 Gauss que es 10 3 veces superior a los campos 
que se deducen que tengan los radio pulsares. Entonces, si el efecto dinamo opera globalmente 
dependera de la relation rotacion-conveccion. Asi una ENs con perfodos de rotation del orden o 
superiores al perfodo de conveccion de 10 milisegundos lograra campos de 10 15 Gauss. Ellas son 
las llamadas Magnetars (Magnetic Stars-Estrellas Magneticas). Por otro lado el Efecto Dinamo 
no tendra lugar en pulsares como el Pulsar del Cangrejo que rota una vez cada 20 milisegundos, 
mucho menor que el perfodo de conveccion. En ese caso el campo magnetico no se incrementarfa 
mas alia de 10 12 Gauss. 

Hemos visto dos mecanismos que de conjunto podrian explicar campos magneticos intensos 
en las ENs llamadas Magnetars. No obstante a pesar de que estas teorfas son reconocidas por la 
comunidad cientifica no se descartan otras ideas para explicar los colosales campos magneticos 
que en ENs aparecen [21] . 

1.3. Estrellas de Quarks. 

El primer modelo de EQs fue propuesto por Itoh en el ano 1970 [22], con este modelo simple 
se investigo el equilibrio hidrostatico de las EQs siguiendo el metodo de Oppenheirmer y Volkoff. 
La MEQ puede encontrarse como una fase en el nucleo de las ENs en cuyo caso tendremos una 
estrella hfbrida, o puede formar una EQs, es decir, una estrella formada en su totalidad por 
quarks con una pequena capa exterior de electrones. 
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Debido a la presencia de la fuerza nuclear fuerte las EQs son sistemas autoligados, a di- 
ferencia de la materia nuclear, que esta ligada mediante la gravedad. Esto hace que las EQs 
puedan presentar configuraciones mas compactas, es decir, menores radios con igual masa. 

Las EQs al ser autoligadas no presentan una masa minima como las ENs (M m i n = 0,1 M ) [3J , 
y su masa maxima tiene una dependencia con el radio maximo (para masas pequenas en com- 
paracion con la masa maxima) del tipo M ~ R 3 , en contraste con las ENs, cuyo radio decrece 
con el aumento de la masa (M ~ R~ 3 ). 

Otra caracteristica que distingue a las EQs es que violarian el lfmite de Eddington [UJ. El 
lfmite de Eddington implica la existencia de una luminosidad maxima para una estrella con una 
masa dada, tal que la configuration se encuentre en equilibrio hidrostatico. Para una estrella 
normal, el lfmite de Eddington esta dado por 

L Edd ~ 8,1 x 10 43 ^- MeVs" 1 . 

M Q 

Para cualquier valor de radiation que supere este lfmite, no habra equilibrio hidrostatico, 
causando la perdida de masa de la estrella. El mecanismo de emision en una EQs producirfa 
luminosidades por encima de dicho lfmite. Al ser autoligada la EQs, el lfmite de Eddington 
no le es aplicable. Asf, su superficie puede alcanzar temperaturas mucho mayores que las que 
se observan en una ENs, y por lo tanto, la emision de energfa termica |23j tambien debe ser 
mayor. 

A pesar de estas diferencias experimentalmente es dificil diferenciar a una ENs de una EQs. 



1.4. Evidencias observacionales. 

La primera ENs fue descubierta en 1967 por Jocelyn Bell al estudiar las senales de radio 
provenientes de un objeto en el cielo. Jocelyn descubrio el pulsar CP 1919 (Cambridge Pulsar 
19 h 19 m ) observando que la senal de radio que recibfa tenia un perfodo exacto de 1.3373 segundos 
y un ancho de pulso de 0.04 segundos [21]. Al inicio se penso incluso que esta senal provenfa 
de alguna civilization extraterrestre, hasta que Thomas Gold y Fred Hoyle identificaron que la 
senal se generaba a partir de una ENs en rotacion con un elevado campo magnetico [25] . 

En la actualidad han sido identificados mas de 2000 objetos como ENs. La determinacion de 
parametros como la masa, el radio y las caracterfsticas de enfriamiento pueden imponer restric- 
ciones a las EdE para la materia a altas densidades. Los astronomos buscan fundamentalmente 
estrellas aisladas para determinar sus propiedades, a partir de los espectros de emision, sin que 
interfieran complicados mecanismos de transmision de materia de un objeto a otro [26J. 

La mejor forma de determinar la masa de una ENs es a traves de la tercera ley de Kepler [3]. 
Uno de los mejores resultados en la determinacion de masas se obtuvo para el pulsar binario 
PSR 1913+16 cuyas masas son 1,3867 ± 0,0002 y 1,4414 ± 0,0002 M [27] respectivamente. 

La determinacion de los radios se hace de forma indirecta, se combinan los datos observa- 
cionales con la teorfa para obtener algunos estimados que arrojan valores entre 9 — 15 km [llj. 

El campo magnetico se determina a partir de las observaciones del decrecimiento en el 
perfodo de rotacion de los pulsares. Aunque estas mediciones no son muy precisas, se estiman 
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valores de entre 10 12 y 10 15 G (28). 

El corrimiento al rojo z s se determina a partir del analisis del espectro de emision de la 
estrella. 

Existe toda una serie de observaciones que hacen pensar que algunos de los objetos que se 
comportan como ENs, son en realidad EQs o estrellas hibridas. 

El observatorio Chandra de rayos X de la NASA encontro dos estrellas inusuales: la fuente 
RX J1856. 5-3754 [29l 130] con una temperatura de 10 5 K y la fuente 3C58 con un perfodo de 
65 ms. RX J1856. 5-3754 es demasiado pequena para ser una ENs convencional y 3C58 parece 
haberse enfriado demasiado rapido en el tiempo de vida que se le estima. 

Combinando los datos del Chandra y del telescopio espacial Hubble, los astronomos deter- 
minaron que RX J1856. 5-3754 radia como si fuera un cuerpo solido con una temperatura de 
unos 7 x 10 5 °C y que tiene un diametro de alrededor de 11 km, que es un tamano demasiado 
pequeno como para conciliarlo con los modelos conocidos de las ENs. 

La fuente de rayos X denominada EXO 0748-676 es otro de los objetos que atrae la atencion 
pues posee una masa de 1,5 < M/M < 2,3 y un radio de 9,5 km < R < 15 km con un 
corrimiento al rojo de z s = 0,35 |31| . 

En febrero de 1987 en el observatorio de Las Campanas, en Chile, fue observada la Supernova 
1987A en la Gran Nube de Magallanes. El remanente estelar que ha quedado como consecuencia 
de la explosion de esta Supernova, podrfa ser una EQs, ya que el periodo de emision de este 
pulsar es de P em = 0, 5ms. Una ENs canonica no podrfa tener una frecuencia de rotation tan 
alta P21E2]. 
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2.1. Modelo para describir la materia de quarks. 

2.1.1. CDC y materia de quarks. 

El Lagrangiano de la CDC es |33j: 

Cqcd = q(il»D,-m)q - G\ v , (2.1) 

donde q es el campo de los quarks con seis sabores (u,d,s,c,b,t) y tres colores, y m = 
diagf(m M , m^, . . . ) es la matriz de masa. La derivada covariante 



D, = d„ - ig ^ Al (2.2) 



esta relacionada con el campo de los gluones A a ^ y 



G%, = d,A a u - d v A; + gf^A^Al (2.3) 

es el tensor de fuerzas de los gluones. A a y f abc son los generadores del grupo SU(3) (matrices 
de Gell-Mann) y las correspondientes constantes antisimetricas de estructura respectivamente. 
g es la constante de acoplamiento de la CDC. 



El Lagrangiano (2.1) es por construccion simetrico bajo las transformaciones del grupo 
SU (3) en el espacio de color. Este grupo tiene un caracter no Abeliano debido a la presencia de 
los terminos f abc en (2.3). Debido a esto, como ya se explico en la introduccion, la CDC presenta 



algunas caracterfsticas peculiares, que no estan presentes en teorfas de calibracion Abelianas 
como la Electrodinamica Cuantica, ellas son: 

■ Los gluones presentan carga de color. 



10 
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La CDC es una teoria asintoticamente libre, es decir el acoplamiento disminuye para 
cortas distancias o de forma equivalente para grandes energfas. En la aproximacion de un 
lazo tenemos que: 

= ^ 4?F (9 A) 

" s = 4^ (ll-lN f )ln(^/A* CDc y { ■ } 

donde Nf es el numero de sabores y Acdc — (200 ~ 300) MeV es un parametro de 
renormalizacion. 

Para bajas energias la interaction es fuerte, es decir, la fuerza entre las particulas aumenta 
con la distancia (confmamiento). 

El Lagrangiano presenta (aproximadamente) simetrfa quiral, es decir, es simetrico ante 
las transformaciones globales SU(Nf)L x SU(Nf)n. Esta simetrfa serfa exacta en el lfmite 
en que el numero de sabores de quarks Nf, tenga masa nula. 



Veamos con mas detenimiento las implicaciones de la ecuacion (2.4) [34]. En el caso de 
la Electrodinamica la constante de acoplamiento (constante de estructura fina) tiene un valor 
olem — 1/137 < 0,01, este hecho hace que podamos dar un tratamiento perturbativo a esta 
teoria. 

Para la CDC se podrian utilizar metodos perturbativos en el limite /i — > oo, esto equivale, 
para la MEQ a la condition de que la densidad barionica tienda a infinito (% — >■ oo) ya que el 
potencial qumiico barionico seria fx ~ 

Estimemos el valor de a s para las densidades que se encuentran en el interior de las EQs. 
Para densidades tib = lOrio, obtenemos que a s ~ 0,8. Incluso si ns = 10 6 rio obtenemos que 
a s ~ 0,15. Para que a s ~ cxem se necesitaria que tib > 10 123 no- 

Estas estimaciones nos aseguran la imposibilidad de utilizar Teoria de Perturbaciones a 
partir de la CDC si intentamos describir un sistema real, como es el interior de los OC. Por 
ello para estudiar estos sistemas se han propuesto varios modelos fenomenologicos, uno de los 
mas exitosos es el Modelo de Bag del MIT, que describiremos seguidamente y que sera el que 
usaremos en esta tesis. 



2.1.2. Modelo de Bag del MIT. 

El modelo de Bag (bolsa) del MIT fue propuesto en los anos 70 [35] para describir desde 
el punto de vista microscopico a los hadrones. En este modelo se logra reproducir la libertad 
asintotica y el confmamiento de los quarks, dos de las propiedades mas importantes de la CDC, 
a traves de un parametro libre del modelo Bh ag . 

En este modelo se describe a los quarks como particulas libres (o casi libres) confinados 
dentro de un Bag. La estabilidad del Bag se garantiza introduciendo el parametro Bb ag que 
se interpreta como una contribucion positiva a la densidad de energfa, y negativa a la presion 
dentro del volumen. De forma equivalente podemos atribuir el termino —Bb ag a la region fuera 
del Bag, lo que conduce a un vacio con una densidad de energia negativa e vac = —Bf, ag y una 
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presion positiva p vac = +Bb ag - El tensor energia impulso dentro del Bag para un fluido perfecto 
viene dado por: 

T" v = diag(e,P,P,P), (2.5) 

donde e es la densidad de energia y P es la presion del fluido. 

El termino B\, ag g^ u (siendo g^ v el tensor metrico) se anade al tensor energia impulso dentro 
del Bag quedando entonces las siguientes expresiones para la presion y para la densidad de 
energia: 

/ 

e = Y^e f + B bag , (2.7) 
/ 

donde la suma es sobre los sabores de los quarks. 

El modelo fenomenologico de Bag del MIT es ampliamente usado para describir la materia 
en el interior de los OC. Es el que utilizaremos en el Capitulo [3] para encontrar las EdE de la 
MEQM y estudiar la estabilidad de la misma. 



2.2. Hipotesis de la materia extrana de Quarks. 

La hipotesis sobre la existencia de la MEQ o conjetura de Bodmer-Witten establece que la 
materia formada por quarks u, d y s en estado libre, puede ser el estado basico de la materia. 
Este hecho se cumple si se satisface la desigualdad: 



E 
A 



E 

< J 

MEQ 



> 

56 Fe 



(2- 



es decir, la energia por barion para la MEQ es menor que la energia por barion para el 56 Fe 
(nucleo mas estable de la naturaleza). 

El hecho empfrico de que en la Naturaleza la materia que se observa esta formada por nucleos 
compuestos por neutrones y protones, y que no hallemos la MEQ, no entra en contradiccion con 
la hipotesis de su existencia. Para que esta se forme a partir de nucleos normales es necesario 
que los quarks ujd, que forman los neutrones y protones, se transformen en quarks s. Pero no 
solo con eso basta, se necesita para lograr la estabilidad que la fraccion de quarks s sea grande 
(n s ~ n u ~ rid). La probabilidad de que ocurra la reaccion debil que garantiza que se formen 
los quarks s es casi nula, lo que hace que solo para grandes densidades, en el interior de las 
ENs y durante millones de anos, se forme la cantidad necesaria de quarks s para que la MEQ 
sea mas estable que la materia nuclear [36]. 

Sin embargo, la estabilidad de los nucleos si excluye la estabilidad de la materia compuesta 
solamente por quarks u y d, como se demuestra a continuacion. 



Capitulo 2: Materia de quarks. 



13 



2.2.1. Estabilidad de la materia de quarks normal en equilibrio este- 
lar. 

Analicemos la estabilidad para el caso de la materia de quarks formada por quarks u y d en 
equilibrio /3 con electrones. Este tipo de materia pudiera formarse en el interior de los OC al 
formarse el plasma quark-gluon. Como veremos, la misma constituirfa un estado metaestable, 
pues tiene mayor energfa por barion que la materia nuclear. 

El equilibrio qufmico dentro de la estrella se garantiza a traves de la reaction de desintegra- 
cion j3: 

d <-> u + e~ + v e (2.9) 

Como las estrellas son neutras se impone la condicion de neutralidad de carga que en el caso 
general se expresa mediante: 

J2n f q f -n e = 0, (2.10) 

/ 

donde qf (/ = u,d,s) son las cargas de los quarks que vienen dadas en la Tabla[T]y (i = 
u, d, s, e) (emplearemos el submdice / cuando tengamos en cuenta solamente los quarks y el 
submdice i cuando incluyamos a los electrones) son las densidades de las partfculas involucradas. 
Se impone tambien la conservacion del numero barionico: 

n B = (2.11) 

De las condiciones anteriores podemos escribir el siguiente sistema ecuaciones: 

Hu + ^e — ^d = equilibrio (3 (2.12a) 
2n u — rid — 3n e = neutralidad de carga (2.12b) 
n u + n<i — 3riB = conservacion del numero barionico, (2.12c) 

donde fi^ y n; son los potenciales qumiicos de las partfculas. 

Hemos supuesto que ha pasado el tiempo suficiente para que los neutrinos hayan escapado 



del sistema por lo que tomaremos fi Pe = 0. El sistema (2.12) es un sistema de tres ecuaciones 
con tres incognitas fii (i = u,d,e). Los resultados obtenidos a partir de este sistema seran 
utilizados para comparar con los de la MEQ. 

Podemos realizar un analisis de la estabilidad de este tipo de materia. Para esto tendremos 
en cuenta que en equilibrio /3 la fraccion de electrones es despreciable, tambien consideraremos 
que las masas de los quarks u, d son nulas. En este caso tenemos: 

,3 Q.,4 ,,4 



7T 2 / 47T 2 f 47T 2 

Utilizando el modelo de Bag del MIT e imponiendo la condicion de equilibrio (P = 0) 



obtenemos para la energfa y la presion total dentro del bag (2.6) y (2.7) 
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/ 

e = 4B bag , (2.15) 

La neutralidad de carga en este caso toma la forma: rid = 2n u =>- /i„ = 2~ 1 / 3 fid = yU^, 
para la presion tenemos: P u d = P u + Pd = (^f^Vud = y la conservacion del numero 
barionico nos da: n Bu d = riu 1 j nd = fr u d/^ 2 , con estas ecuaciones podemos calcular la energfa por 
barion para el gas formado por quarks u y d: 

= (1 + 2 4 / 3 ) 1 / 4 (4 7 r 2 ) 1 / 4 J B fe 1 4 4 = 6,4115^ ~ 943 MeV (2.16) 

donde hemos tornado B l J a A g = 145 MeV =^ B bag ~ 57 MeV fm -3 . La energfa por barion 
para un gas de neutrones es la masa del neutron (m n = 939, 6 MeV), para un gas de 56 Fe 
la energfa por barion se puede calcular como: E/A\ha-pe — (56mAr — 56 ■ 8,8)/56 = 930 MeV. 
Por tanto la estabilidad del gas formado por quarks u y d con respecto al gas de neutrones 
requiere que tudl^Bud < m n =>• Bbag < 60MeVfm~ 3 y con respecto al gas de 56 Fe que 
e ud /n Bud < E/A\ 5&Fc B bag < 57MeVfm" 3 . 

En la naturaleza se observan los neutrones y el 56 Fe pero no el gas de quarks u y d, por lo 
que se concluye que el B bag debe ser mayor que los valores antes mencionados. Este hecho hace 
que se tome el valor B bag = 57MeVfm~ 3 como un lfmite inferior para este parametro. 

2.2.2. Estabilidad de la Materia Extrana de Quarks en equilibrio es- 
telar. 

Analicemos ahora la estabilidad de la MEQ, es decir, materia formada por quarks u, d y s 
libres. El equilibrio qufmico de este sistema esta garantizado por las interacciones debiles: 

d <r> u + e _ + u e 
s -H- u + e~ + z/ e 
U + S -H- d + u 

estas reacciones conducen a relaciones entre los potenciales qufmicos. Tomaremos igualmente 
que /i Pe = pues los neutrinos escapan de la estrella. El equilibrio (3, la neutralidad de carga y 
la conservacion del numero barionico quedan expresados mediante: 

H u + yL e - n d = , fid-l^s = 0, (2.17a) 
2n u - n d - n s - 3n e = 0, (2.17b) 
n u + n d + n s - 3n B = 0. (2.17c) 

Al igual que en el caso de la materia de quarks normal supondremos que la fraccion de 
electrones es muy pequena, por tanto la neutralidad de carga queda como: n u = rid = n s =>• 
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l^uds = = = Us- Para un B^ ag fijo el gas de u , d , s ejerce la misma presion que el gas de u y 
d, tambien la energia es la misma (e u d s = ABi, ag ). Esto ocurre cuando fx u d s = 3~ 1//4 (l+2 4/ ' 3 ) 1//4 /x 11( 2 
con lo que se obtiene el numero barionico: riBuds — Huds/ft 2 = [(1 + 2 4//3 )/3] 3//4 n j B U d ; ssto nos da 
una energia por barion: 

= 3 3 / 4 (4 7 r 2 ) 1 / 4 J B 6 1 a / o 4 = S^lBl 1 ^ ~ 829 MeV, (2.18) 

n-Buds 

este valor es 100 veces menor que para la materia normal de quarks (para un mismo valor de 
Bbag)i por lo que podemos decir que se cumplen las siguientes desigualdades: 



E 




E 






< 






< a 

MEQ ^ 


56 Fe 


A 


u 



(2.19) 



donde el parametro Bb ag se encuentra entre 57 — 90 MeV fm~ 3 . Este resultado conduce 
directamente a la hipotesis de Bodmer-Witten para el caso de ambientes ricos en quarks s 
como puede ser el interior de las ENs [37J. 

En el capitulo siguiente estudiaremos los efectos del campo magnetico en la estabilidad de 
la MEQ y de la materia normal de quarks. 



Capitulo 3 



Propiedades termodinamicas de la 
Materia Extrana de Quarks. 



3.1. Materia extrana de quarks a5 = 0yT = 0. 

De forma general el gran potencial termodinamico para el gas de quarks relativistas y libres, 
en el marco del modelo de Bag, viene dado por cuatro terminos 



n = B bag V + £ [Q hV (fi t ,T)V + Qi, s (^,T)S + QiA^,T)C] , (3.1) 

i 

donde la suma se extiende sobre i = (u, d, s, e), y tiene en cuenta todas las especies de fermiones 
presentes: electrones, quarks u, d y s; fa es el potencial qufmico de cada gas de partfculas; T la 
temperatura; V, S y C son el volumen dentro del cual se encuentra el gas, el area y curvatura 
de la superficie del Bag respectivamente. 

Consideraremos en nuestro estudio un volumen en el interior de la estrella por lo que no 
tendremos en cuenta los efectos superficiales ni de curvatura. Dentro de este volumen, los quarks 



u, d, s se comportan como particulas libres Fig (3.1 ). El termino volumetrico del Gran Potencial 



Termodinamico para fermiones libres esta dado por la expresion 

Qi,MiT) = jMf+{p,to,T)rfalk,T))<Pp], (3.2) 

T) = l + exp{-^^}, (3.3) 

aquf / representa las contribuciones de las partfculas y antiparticulas (/ + para particulas y f~ 
para antiparticulas), E p ^ es la energfa de las partfculas dada por el espectro y di es el factor de 
degeneracion: 1 para electrones y 3 para quarks^] El Gran Potencial juega un papel importante 



toma por convenio las unidades naturales h ■ 
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Figura 3.1: Volumen local en el cual los quarks se comportan como particulas libres. 



ya que de el se derivan todas las propiedades termodinamicas del sistema ED] y se pueden 
obtener las EdE del sistema. 



En (3.2) y (3.3), E p ^ se obtiene de la solution de la ecuacion de Dirac para fermiones libres: 



E Pti = ^P 2 + ml (3.4) 

donde pyrrii representan el momentum lineal y la masa en reposo de cada partlcula respecti- 
vamente. 

En nuestro caso podemos hacer un cambio de variables a coordenadas esfericas e integrar 
por partes con lo que obtenemos: 

a>„T) = -A r / ^±_ (n F (^^) + np(^±^)), (3.5) 
67r Jo \/p 2 + m 2 V T T J 

donde uf(x) = l/(e x + 1) es la funcion de distribucion de Fermi-Dirac. 

En el caso de una ENs o una EQs la temperatura es mucho menor que la temperatura de 
Fermi (T/Tp ~ 10 -4 ), por lo que podemos considerar el lmiite degenerado (T = 0). Esto hace 
que np = 8(fii — E p j) donde 9{x) es la funcion paso unitario, ademas la contribucion de las 

antiparticulas es nula pues np( Ep 'ljt^ ) > con lo que obtenemos: 

<M/*.0) = -A / (3.6) 

o v 
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La integral (3.6) da como resultado: 



5 3 

toi,v(jH,0) = -^oixiix^ - -)p F ,i + -cosh" 1 ^]}, 



(3.7) 



d 'TTb- 

donde hemos definido la variable adimensional Xi = fii/rrii, Sl — ^ \ J VF,i — a/x? — 1 es el 

momentum de Fermi adimensional de cada especie. 

A partir de (3.7) podemos obtener la presion, la densidad de partfculas y la energfa del 
sistema: 



P 

Ni 



i 

<9a 



(3.8) 

(3.9) 
(3.10) 



la energfa total viene dada por la ecuacion (2.7). Se recuperan las expresiones (2.13) cuando 
tomamos m,- = 0. 



3.2. Materia Extrana de Quarks Magnet izada a T = 0. 

En el caso de la MEQM, estudiaremos igualmente un volumen local con un campo magne- 



tico, que consideraremos constante, homogeneo y en la direccion del eje z Fig (3.2). 

El espectro energetico se obtiene al resolver la ecuacion de Dirac para fermiones en presencia 
de un campo magnetico y tiene la forma |41j: 




E n J = \ P 2 Z +P\+ ml p 2 ± = mf{( J—(2n + 1 - V ) + 1 - ^B) 1 - 1}. (3.11) 



qi mi 2n 2m e 

donde n son los niveles discretos de Landau, 77 = ±1 son los autovalores del operador de spin, 
B es el campo magnetico, B1 es el valor del campo critico (campo para el cual la energfa 
ciclotronica se iguala a la energfa en reposo), g^, Qi y representan la carga electrica, el 
momento magnetico anomalo (MMA) y la masa de cada partfcula respectivamente. La carga 
electrica de cada partfcula esta en la Tabla [TJ los valores de MMA y la masa que utilizaremos 
se dan en la Tabla 13.21 

Analogamente al proceso que seguimos para calcular el potencial termodinamico para el gas 



de quarks sin campo magnetico, en el caso de campo tenemos que calcular la integral (3.2). 

En presencia del campo magnetico la simetrfa deja de ser esferica, por lo que para integrar 
realizamos un cambio a coordenadas cilfndricas obteniendo: 



ViA^T) = £ £ / [ln[l + exp{-^-^}]rfpJ, (3.12) 
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Figura 3.2: Igual que en la Fig (3.1), pero considerando un campo magnetico constante y 
homogeneo. 



Particulas 


u 


d 


s 


e 




1,85/i^ 


— 0,97/iAr 


— 0,58/iAr 


1,16 x 10" 3 /ii 


rrii (MeV) 


5 


5 


150 


0.5 



^Magneton Nuclear fx N = e/2m p ~ 3,15 x 10~ 18 MeV G _1 . 
tMagneton de Bohr [i B = e/2m e ~ 5,79 x 10~ 15 MeV G _1 . 

Tabla 3.1: Valores del MMA y masa de las particulas consideradas en nuestro estudio. 



oo p+oo 



oo J — oo 



[dp x dp y ] = 2irqiB. 
El potencial termodinamico tiene la forma: 

,ij2, X i+PF,, 



donde 



i n t/=± 1 



hi 



47r 2 ' / '/'.' ~ \/' r > ~ ^ 2 ' 



B 



B' 



- (2n + 1 - 77) + 1 - r]y t B. 



(3.13) 



(3.14) 



(3.15) 
(3.16) 
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En la ecuacion (3.14) la suma sobre los niveles de Landau n es hasta el valor entero n l max = 
I [{{xi + rn/iB) 2 — 1) Bf/ (2B)], ppi y son el momento de Fermi y la masa magnetica respec- 
tivamente. 

La presencia del campo magnetico hace que surja una anisotropfa en la presion, de tal forma 
que la podamos dividir en paralela P\\ y perpendicular al campo P±, por lo que el tensor energfa 
impulso (2.5) queda: 

T» u = diag(e,P ± ,P ± ,P {l ), (3.17) 

donde Pn = —Q y P± = —Q — M.B siendo Ai = —dVt i v /dB la magnetization. 
Los parametros termodinamicos se expresan mediante las expresiones: 

i 1 n 7]=±1 

^max / I T) \ 

e = Q + vN = Bj2M°J2 E UiP^ + ^ln^pJ, (3.19) 

2 X i + p Fi 



"•max / T ■ A- T) \ 

P\\ = -n = BY J M°J2J: Uph-hf^-^Y (3.20) 

i n ?7=±1 ^ % / 

P ± = -n-MB = Bj2M°Y: E 2^1n^ , (3.21) 

n rj=±l N re * / 



donde 



JV? = (3.22) 

B (2n + 1 - 7?) , 

ll = ; — TjyjB, 3.23 

25,V(2n + 1 - V )B/Bf + 1 

son magnitudes adimensionales. 

Las expresiones anteriores contienen las contribuciones del diamagnetismo de Landau (cuan- 
tizacion de los niveles de Landau), y el paramagnetismo de Pauli debido a la presencia del MM A. 
El paramagnetismo de Pauli no modifica las EdE ni la estabilidad de la MEQM, por lo que no 
tendremos en cuenta su efecto en lo que sigue. 



3.2.1. Ecuaciones de estado para la MEQ. 

Para determinar las EdE de la MEQM en equilibrio estelar, hay que resolver numericamente 
las ecuaciones (2.17). En la Fig (3.3) mostramos las EdE para la MEQ (B = 0) y para la 
MEQM para B = 5 x 10 18 G, en todos los casos hemos tornado el valor Bt, ag = 75MeVfm -3 . 
Para comparar hemos incluido las EdE de la materia normal (u, d en equilibrio con electrones). 
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El caso B = no muestra diferencias con la ecuacion de estado P = (e — 4Bb ag )/3 de la 
materia formada por quarks u y d sin masa, analizada en el epigrafe (2.2.1), lo que muestra 
que considerar las masas de los quark u y d diferentes de cero, siendo ellas tan pequenas, no 
influye de manera significativa en los resultados. 

Aunque las variaciones que el campo magnetico produce en las EdE no son muy notables, 
ellas pueden alterar las propiedades macroscopicas como son la masa y el radio de la estrella 
como veremos en el siguiente capitulo. 




Figura 3.3: Ecuacion de estado para la MEQ: presion contra energia para diferentes valores de 
campo. Para comparar se han anadido las graficas para la materia normal en equilibrio con 
electrones. 



3.3. Estabilidad de la materia de quarks magnet izada. 

Estudiemos la estabilidad de la MEQM como funcion de los parametros del modelo, es decir, 
la densidad barionica (tib), el campo magnetico, el parametro de Bag y la masa del quark s 
m s . Debido a la presencia del campo magnetico, la anisotropia de las presiones implica que 
P± < P|| [13J. La condicion de estabilidad se expresa mediante la relacion: 

p± = J2 p ±,i- b ^ = °- ( 3 - 24 ) 



Para investigar la estabilidad de la MEQM en equilibrio estelar resolveremos el sistema de 
ecuaciones (2.17) y (3.24), utilizando los resultados obtenidos en (3.21)- (3.18). Al resolver 
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este sistema de ecuaciones se obtiene la region en la cual los parametros cumplen con las 
desigualdades de estabilidad: 
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(3.25) 




Figura 3.4: Energia por barion en funcion de la 
densidad barionica para B = y B = 5 x 10 18 G, 
asumiendo la condicion de estabilidad de presion 
cero. 



Figura 3.5: Energia por barion contra presion 
para B = 0y5 = 5x 10 18 G. Tomando B h&g = 
75MeVfm~ 3 . La linea horizontal de puntos co- 
rresponde a E/A\ 5 e Fe ~ 930 MeV. 



En la Fig (3.4) presentamos una comparacion de la energia por barion E/A (equivalente- 
mente s/ub) contra la densidad de partfculas n^/rzo para B = 0y5 = 5x 10 18 G exigiendo 
la condicion de equilibrio P± = 0. Hemos asumido que m u = = 5 MeV y m s = 150 MeV en 
todos los casos. 

Como puede observarse en la figura E/A es menor cuando esta presente el campo magnetico. 
Para B = 5 x 10 18 G obtenemos que E/A « 919 MeV, n B /n ~ 2,2 y B bag ~ 75MeVfm" 3 
mientras que para B = e igual valor del Bag, E/A ~ 929 MeV, ns/n ~ 2,1. 

El comportamiento de E/A con la presion se muestra en la Fig (3.5) fijando Bb ag = 
75MeVfm -3 . Podemos notar que el punto de presion cero para la MEQM se alcanza para 
un valor de densidad de energia menor que en el caso de la MEQ. Por tanto la MEQ es mas 
estable y mas compacta en presencia de campo magnetico. 

Un estudio mas detallado de la estabilidad lo podemos realizar a traves de las ventanas 
de estabilidad, las cuales nos brindan una mayor informacion acerca de los lmiites para los 



parametros entre los cuales se cumple la desigualdad de estabilidad (3.25). Para investigar 



como el campo magnetico afecta estas ventanas, consideraremos las regiones de estabilidad en 
el piano (m s ,n^). Para ilustrar, fijaremos el campo magnetico en el valor B = 5 x 10 18 G, 
B = G y estudiaremos los contornos Bb ag = const y E/A = const. 
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Figura 3.6: Ventanas de estabilidad para la 
MEQ en el piano (m s , ng). Se muestran los con- 
tornos de B^ag = const y E/A = const. 



Figura 3.7: Ventanas de estabilidad para la 
MEQM en el piano (m s , Ub)- Se muestran los 
contornos de Bbag = const y E/A = const. 



Estos resultados se presentan en las Fig (3.6) y Fig (3.7) para MEQ y MEQM respecti- 
vamente. Como puede observarse, el campo magnetico hace que las ventanas de estabilidad 
para la MEQ se muevan hacia valores mayores de la densidad barionica. Ademas se modifica el 
intervalo permitido de valores para el parametro de Bag debido a que por debajo del contorno 
de energia de 930 MeV las EdE para MEQM corresponden a un valor de E/A a presion cero 
menor que el del 56 Fe 



A I MEQ 



< 



E 
A 



56 Fe J 



Del analisis de la Fig (3.7) podemos considerar la MEQM como absolutamente estable si 
3 < Bb ag < 90MeVfm~ 3 , 1,85 < ub/uq < 2,6 y para un amplio rango de masas 

75 MeV fin" 3 



57MeVfm- d < B bag < 90MeVfm 
del quark s, en particular para m s = 150 MeV y Bb ag 
en el proximo capitulo para calcular algunos parametros de las EQs. 



valores que utilizaremos 



Capitulo 4 

Estrellas de Quarks magnetizadas. 



4.1. Observables. 

4.1.1. Relacion Masa Radio y ecuaciones TOV. 

En el capitulo anterior estudiamos la estabilidad y las EdE de un gas magnetizado de quarks. 
^Como podemos vincular esto con la curvatura en el espacio tiempo, causada por las enormes 
densidades de la materia en una EQs? 

El problema se puede tratar separando la fuerza gravitacional que tiene un rango de accion 
infinito y las fuerzas fuertes de corto alcance. De esta manera en nuestro estudio anterior 
hemos considerado que en las escalas de las interacciones fuertes (~ 1 fm) la curvatura es nula, 
la gravitation no es necesaria incluirla. 

Ahora bien, para las escalas de las EQs, (~ 10 km) no podemos obviar el papel de la 
curvatura, es decir, la gravedad curva el espacio tiempo a escalas macroscopicas |12] . 

Por esto para conocer la estructura (radio, masa) de una estrella es necesario utilizar la 
Teorfa General de la Relatividad de Einstein. 

Las ecuaciones de Einstein establecen un estrecho vinculo entre el contenido de materia en 
un lugar del espacio y la curvatura del mismo p]: 

G\ = kT^ (4.1) 

{H,v = 0,1,2,3), k = 8ttG n G n = 1,32 x 10~ 42 fmMeV" 1 , G M „ = - \Rg^ u es el tensor de 
Einstein el cual viene determinado por el tensor de Ricci R^ v y por el escalar de Ricci R = R 11 ^, 
estos ultimos dependen de segundas derivadas de la metrica. 

p pa pa I pa p/3 p/3 pa f a n\ 

^ixv — 1 1 iia,u 1 fiu L a p L fia 1 v$i K^-^J 

las cantidades T° u son los indices de Christoffel, que dependen de primeras derivadas de la 
metrica por la formula, 

T ^u = —(gpw + 9vp,n - 9nv,p)- (4.3) 



2 



El tensor energia impulso T 7 ^ viene dado por la expresion (3.17). 
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Para una estrella estatica en equilibrio se emplea la metrica: 

ds 2 = - e ^(r) dt 2 + e A(r) dr 2 + ^2 + r 2 ^2 



(4.4) 



Utilizando la metrica (4.4), la ecuacion (4.1 ) y la ley de conservation de la energia (T^ = 0) 
podemos encontrar las ecuaciones TOV [lj, que nos dan las configuraciones de estrellas estaticas 
y con simetrfa esferica: 



dM 

dr 
dP 

dr 



47rGe(r) 

(e(r) + P(r))(M(r) + AnP(r)r 3 ) 



-G- 



2rM(r) 



(4.5a) 
(4.5b) 



Para resolver el sistema de ecuaciones (4.5) utilizaremos las soluciones numericas de las 



EdE para la MEQM (3.18) y (3.21). El radio R y la masa correspondiente M de la estrella 
se determinan imponiendo la condition de presion cero P(R) = 0. La presion central queda 
fijada por la EdE, -P(O) = P c , si imponemos ademas la condition M(0) = podemos resolver 
el sistema de ecuaciones TOV. 



i 
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i 
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B = 0G 
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n B /n Q 



Figura 4.1: Configuraciones de M — R obtenidas 
para MEQM. Se muestran curvas para diferentes 
valores de campo magnetico B. Para comparar 
se han anadido las graficas M — R para la ma- 
teria de quarks normal. Las configuraciones mas 
compactas son las de MEQM 



Figura 4.2: Masa contra densidad central para 
diferentes valores del campo magnetico. Se mues- 
tran los resultados para el gas formado por u y 
d. Se puede observar que el aumento del cam- 
po tiende a hacer decrecer la masa maxima para 
una densidad fija. 



En la Fig (4.1) mostramos los graficos de las configuraciones de equilibrio para la MEQ 
y para la MEQM con un valor del parametro de Bag {B^ ag = 75MeVfm -3 ). Como ya fue 
mencionado en el capitulo anterior, el campo magnetico hace que la MEQ sea mas estable por 
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lo que las estrellas pueden ser mas compactas. Se observa como a medida que aumenta el campo 
las configuraciones obtenidas presentan menores radios. 



La solution de las ecuaciones TOV (4.5) representadas en la Fig (4.1 ) dan las configuraciones 
de equilibrio para la estrella, pero este equilibrio puede ser estable o inestable. 

Para conocer si una configuration es estable o inestable hay que estudiar la dependencia 



de la masa con densidad central de masa p c (equivalentemente %) Fig (4.2). La estabilidad 
de la estrella puede determinarse por la pendiente de esta curva. Cuando dM/dp c < las 
configuraciones de equilibrio se consideran inestables [42J. 



En la Tabla 4.1 se muestran los valores de masa maxima con su correspondiente radio y 



densidad central para diferentes valores de campo obtenidos de la Fig (4.1) y Fig (4.2). 



Materia 


M max (M/M ) 


R (km) 


B(G) 




MEQM 


1.56 


8.62 


1 x 10 19 


9.46 


MEQM 


1.62 


8.91 


5 x 10 18 


9.14 


MEQ 


1.65 


9.07 





8.48 


ud 


1.69 


9.20 


1 x 10 19 


7.88 


ud 


1.74 


9.50 


5 x 10 18 


7.31 


ud 


1.77 


9.60 





7.02 



Tabla 4.1: Resultados obtenidos de la Fig (4.1). 



4.1.2. Masa Barionica. 



La masa que encontramos en el epigrafe anterior al resolver las ecuaciones TOV (4.5) es la 
masa gravitacional (Mq)- Hallemos ahora la masa barionica [43j. Ella difiere de la gravitacional 
porque vamos a considerar que la densidad barionica depende de la metrica. 

Para hallarla tengamos en cuenta que el volumen de una capa esferica de la estrella en la 
metrica de Schwarzschild esta dada por la expresion: 

Anr 2 dr 

dV = =5* — — , 4.6 

el numero de bariones dentro de la estrella lo podemos calcular como: 

N B = / n(r)dV = / or - (4.7) 

donde n(r) es el numero de bariones por unidad de volumen. Por tanto la masa barionica de la 
estrella es: 

M B = m n - N B . (4.8) 
Se toma m n como la masa barionica m n = 939 MeV. 
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1.5 



udB = OG 

■ • udB = 5X10 ,8 G 
- -SQM B = OG 
-■-■SQMB = 5X10 1B G 




1.2 



M B /M 

Figura 4.3: Masa barionica contra masa gravi- 
tational. 




Figura 4.4: Masa barionica contra radio. 



La Fig (4.3) muestra la masa gravitatoria en funcion de la masa barionica. La Fig (4.4) 



representan las configuraciones de equilibrio de masa barionica y radios. Como se puede observar 
la masa barionica es siempre superior a la masa gravitacional. La estabilidad puede ser estudiada 
de la misma manera en que se analiza la estabilidad para la masa gravitacional |32]. Podemos 



apoyarnos en la Fig (4.3) que como se puede observar alcanza un maximo y luego comienza a 



decrecer hacia menores valores de masa barionica, esta es la zona de inestabilidad. 



4.1.3. Corrimiento al rojo gravitacional. 

El corrimiento al rojo se define como: 



1 + Z= W> i + z= / ^(receptor) (4g) 

/(fuente) V 9tt(^nte) 

donde gu es la componente temporal de la metrica. 

Esta magnitud da la diferencia de frecuencia de la luz emitida por la superficie de la estrella 
con respecto a la que se mide en un sistema en reposo respecto a la fuente. En este caso estamos 
interesados en el corrimiento al rojo gravitacional, es decir, el cambio de frecuencia que sufre 
la radiacion al moverse en un campo gravitatorio muy intenso como el que existe en las EQs. 



Para el caso de la metrica (4.4) se puede obtener que 

1 



1 - 2G N M/R 



(4.10) 



El corrimiento al rojo se utiliza para descartar EdE que no reproduzcan los datos observados. 
En la Fig (4.5) mostramos los resultados del corrimiento al rojo gravitacional. El valor de 
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Figura 4.5: Corrimiento al rojo gravitacional para distintos valores de campo magnetico y 
Bbog = 75 Me V frrT 3 . 



z s ~ 0,345 obtenido experimentalmente en las lineas espectrales de rayos X proveniente del 
objeto EXO 0748-676 se correspondera en nuestro caso con valores de M/M Q ~ 1,42 resultado 
que esta dentro de los limites de masas medidas experimentalmente. 



Conclusiones. 



En este tesis se "tiende el puente" entre el estudio microscopico de la materia de quarks 
magnetizada (dado por el estudio de las EdE y la estabilidad) y los observables macroscopicos 
de las EQs derivados de las EdE de ese tipo de materia magnetizada. 

El estudio se hace tanto para la MEQ, como para la normal en equilibrio estelar en presencia 
de un campo magnetico intenso, utilizando el modelo fenomenologico de Bag del MIT. 

La estabilidad se ha estudiado con el proposito de conocer si el campo magnetico contribuye 
o no a que sea mas estable la materia de quarks. Las EdE se han obtenido para utilizarlas al 
resolver las ecuaciones de equilibrio hidrodinamico TOV y obtener como los observables (masa 
gravitacional, radios y corrimiento al rojo), se modifican con el campo magnetico. 

Los resultados obtenidos pueden resumirse de la siguiente manera: 

1. Se obtuvo numericamente la EdE para la MEQM Fig (3.3) y para la materia normal de 
quarks en equilibrio estelar. 

2. Se estudio la estabilidad de la MEQM en equilibrio estelar: 



Se verified que las desigualdades de estabilidad (3.25) se cumplen para un amplio 
rango de parametros de modelo de Bag. 

La MEQM es mas estable que la MEQ y que la materia de quarks normal. 

Se obtuvieron las ventanas de estabilidad: curvas en el piano n/ns contra m s de la 
MEQM teniendo en cuenta la variation de la masa del quark s, la densidad barionica, 
el campo magnetico y el parametro de Bag. 

Mientras que para la MEQ el rango de la densidad barionica permitido esta entre 
1)8 ^ ns/no < 2,4 para masas del quark s en el rango 50 < m s < 300 MeV, la 
MEQM admite densidades en el rango 1,85 < ub/uq < 2,6 para rangos de la m s 
50 < m s < 300 MeV y campo magnetico de 5 x 10 18 G. El rango del parametro de 
Bages57< B ba <90MeVfm 



-3 



3. Se resolvieron las ecuaciones TOV con las EdE obtenidas. La relation masa-radio (MR) 
muestra como el campo magnetico contribuye a soluciones estables de EQs mas compactas 
con menores valores de masa gravitacional y radio. 

■ Se obtuvo los efectos del campo magnetico en la masa barionica. Se resolvieron las 
ecuaciones TOV para ella y se obtuvieron las relaciones masa barionica en funcion 
del radio. 
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Conclusiones. 



■ Se estudio el corrimiento al rojo de las EQs magnetizadas. 



Recomendaciones . 



Estudios teoricos recientes han mostrado que a densidades muy altas la materia de quarks 
podria estar en una fase superconductora de color: Color Flavor Locked (CFL) siendo este el 
estado fundamental de la materia [44J . Estudiar como influye la presencia del campo magnetico 
en la estabilidad de esta fase resulta un tema relevante que ya hemos emprendido. 

Por otro lado los observables astrofisicos aquf discutidos merecen que sean comparados 
con datos observacionales. Es esta otra direction de trabajo en que pensamos dirigir nuestros 
esfuerzos. Estudiar otros observables como el momento de inercia de la estrella es tambien algo 
que ha quedado pendiente. 

La inclusion de los efectos de la rotation de las EQs es un tema digno de analizarse, los 
efectos del campo magnetico pueden en este caso resultar importantes. 
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